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Bipolaron formation in organic semiconductors at the interface 
with dielectric gates 
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PACS 73 . 20 . -r - Electron states at surfaces and interfaces 
PACS 7 1 . 38 . Mx - Bipolarons 

PACS 72.80.Le - Polymers; organic compounds (including organic semiconductors) 

Abstract -The formation of the electron-phonon induced bipolaron is shown to be feasible in 
organic semiconductors at the interface with dielectric gates due to the coupling of the carriers 
with interface vibrational modes and to the weak to intermediate strength of bulk electron-electron 
interaction. The polaronic bound states are found to be quite robust in a model with realistic 
strengths of electron coupling to both bulk and interface phonons. The crossover to nearly on-site 
bipolarons occurs for coupling values much smaller than those for nearly on-site polarons, but, 
on the other hand, it gives rise to an activated behavior of mobility with much larger activation 
energies. The results are discussed in connection with rubrene field-effect transistors. 



Introduction In the last years, field-effect devices have 
been realized with correlated oxides and organic materials. 
[T] This has allowed to inject charge carriers into complex 
materials inducing new phases and properties. On the 
other hand, the devices work with solid or liquid dielec- 
tric gates and new effects has to be taken into account: 
for example, the lattice mismatch between bulk and gate, 
the dielectric breakdown for strong electric fields, and the 
interaction of bulk electrons with gate degrees of freedom. 

Single crystal organic field-effect transistors (OFET) arc 
important not only for applications, but also for the under- 
standing of the intrinsic properties of organic semiconduc- 
tors (OS). [5] The most promising and studied systems are 
those realized with oligoacenes, such as rubrene. [3] In the 
OFET geometry, the nearby dielectric can have a strong 
influence on the transport properties. Indeed, some ex- 
periments have pointed out the existence of scaling laws 
of the mobility as a function of the dielectric constant of 
solid [H[S] and liquid [5] gates. If the difference between 
the static and high frequency dielectric constant of the 
gate is small (for example for parylene), the mobility \i of 
these systems exhibits a power- law behavior (fi ~ T~ 7 , 
with 7 ~ 2) at temperatures close or higher than lOO-ft'. 
[3] This conduction process cannot be simply ascribed to 
band transport, [TUTU] and it is believed to arise from the 
interaction of charge carriers with bulk modes. [TOlfTTj On 
the other hand, if the dielectric constant mismatch is high 



(for example for Tt^Os), an insulating behavior is found 
much below room temperature with a strong decrease of 
the mobility. 

The explanation of these behaviors has been based on 
the possibility that the injected charge carriers undergo 
a polaronic localization due to the interaction with the 
polarizable gate dielectrics, leading to the formation of 
Frohlich polarons. [SJH^HH] However, the predicted val- 
ues of the activation energy are smaller than those es- 
timated by accurate calculations. [TS] For example, the 
barrier of rubrene transistors with a Ta^O^ oxide gate is 
calculated to be 12.4 meV, a value much lower than the 
energy of 55 meV found in previous estimates. [5] For this 
reason, recently, the interplay between charge carrier cou- 
pling with bulk vibrations and the long-range interaction 
induced at the interface with the dielectric gate has been 
investigated. |16j These interactions stabilize a polaronic 
state with specific properties, that, however, requires val- 
ues of the bulk electron-phonon (el-ph) coupling still larger 
than those typically estimated. One way to solve this in- 
consistency is to include the effect of disorder. [3lfl7] 

Recently, growth and purification techniques of OS are 
ever and ever improving tending to strongly reduce the 
amount of traps in the bulk and at the interfaces. [18] 
Therefore, in this work, we focus on the clean limit propos- 
ing that the bipolaron formation induced by the long- 
range interaction of the carriers with interface vibrational 
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modes can affect the transport properties. The bipolaron 
mechanism has been invoked in organic single molecules 
[T9] and disordered systems. [2012T] However, bound pola- 
ronic states have not carefully investigated in clean single- 
crystal OS. The polaronic bound states are studied in a 
model, which, in addition to bulk and interface el-ph cou- 
plings, includes local Hubbard and non local long-range 
Coulomb interactions treated through exact diagonaliza- 
tions. The bipolaron formation is possible due to the weak 
to intermediate strength of the bulk electron-electron in- 
teraction in OS. 

In spite of the effects of bulk el-ph interactions, the for- 
mation of bound states takes place for realistic strengths of 
the interface el-ph coupling. With increasing the interface 
interaction, we point out that the bound states form with 
different features: large bipolarons with extended lattice 
distortions for weak to intermediate coupling, small bipo- 
larons with essentially on-site distortions for intermediate 
to strong coupling. In particular, the crossover to small 
bipolarons occurs for coupling values much smaller than 
those for the small polaron formation. On the other hand, 
we will show that the conduction process of small bipo- 
larons is characterized by a hopping behavior with activa- 
tion energies much higher than those of the single small 
polaron, in qualitative agreement with data in rubrene 
OFET. 

Model We study the following Hamiltonian model with 
coupling to bulk and interface vibrational modes |16| , local 
Hubbard and non-local long-range Coulomb interactions: 

H = H el + Hf + H e i-s + + H el _L (1) 



In Eq. ([T]), the electronic part H e i is 
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The kinetic energy operator is in the first term of Eq. ([2} , 
in which t is the bare electron hopping between the near- 
est neighbors on a chain with lattice parameter a, c| a and 
Cj c are the charge carrier creation and annihilation oper- 
ators, respectively, relative to the site Ri and spin a. The 
ab-initio estimate for t is around 50— lOOmeT^. [5] Within 
the model, one considers a one-dimensional stack of pla- 
nar conjugated molecules (axis with the highest hopping 
integral) . [10] The local Hubbard interaction with strength 
U is in the second term of Eq. ([2]), in which the density 
operator of spin a electrons is rii jCr = c\ a Ci, (y . OS are not 
strongly correlated, therefore the quantity U is estimated 
to be less than or of the order of half-bandwidth. [20] The 
non local Coulomb interaction is in the third term of Eq. 
([2]), in which e is the electron charge, e is the bulk relative 
dielectric constant, r\ is a cut-off distance for long-range 
interactions on finite lattices, and n, = J2a n i,o- The ref- 
erence energy for the Coulomb interaction is V = e 2 /ea, 
with e between 2 and 3 for OS. 



In Eq. ([T]), Hg^ is Hamiltonian of free bulk modes. 
For rubrene, the mixing of low frequency intramolecu- 
lar modes with intermolecular slightly dispersive lattice 
phonons takes place. [22j Therefore, we consider a single 
average dispersion-less mode, 23,24 implying 
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where Xi and Pi are the i — th molecule displacement trans- 
verse to the chain and the momentum, respectively, m the 
oscillator mass, k the elastic constant and ujb = \fkjrn 
the mode frequency. These bulk modes are characterized 
by small energies (Tilob — 5 — 10 mcV, with fi Planck con- 
stant). [Hl[To] 

In Eq. (TIJ, H e i-B represents the term similar to the Su- 
Schricffcr-Hceger [23] (SSH) interaction for the coupling to 
bulk modes 



(X i+ l - Xi) (c\ a C i+ l, a + h.c)j , (4) 



with a coupling constant. 1 1 Oil 1 1J This coupling induces a 
modulation of the transfer integral between two neighbor 
molecules. In the adiabatic regime for bulk modes (Tllob <C 
t), the dimensionlcss quantity Xb — a 2 /4kt fully provides 
the strength of the bulk el-ph coupling. The typical values 
of As are in the intermediate regime (in this work, we take 
Xb = 0.1 appropriate for rubrene 

In Eq. ([1]), Hj ^ is the Hamiltonian of free interface 
phonons 

Hf> = Uui ^2 a\a q , (5) 

where w/ is the frequency of optical modes, aj and a q 
are creation and annihilation operators, respectively, rel- 
ative to phonons with momentum q. In comparison with 
models in high dimensions, [13] interface modes are as- 
sumed purely one-dimensional propagating along the elec- 
tron chain lying close to the gate. [TB] 

In Eq. (fTJ), H e i_i is the Hamiltonian describing the elec- 
tron long-range coupling to interface vibrational modes 
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where M q is interaction el-ph term 

M ^ = ghuJlnter piqR] Rp 
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with g dimensionlcss coupling constant, L number of lat- 
tice sites, and i?o cut-off length of the order of the lat- 
tice spacing a. The length Rq indicates the distance of 
the chain from the polar gate. In this work, we take 
i?o = 0.5a and huj] nte r = 0.5t. [T21 [13] In order to 
quantify this coupling, we use the dimensionless quantity 
Xj = J2 q M q /2u!it. We stress that, when one compares 
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this coupling with the more realistic two-dimensonal case, 
one finds a similar long-range behavior. |16j 

In the following, we will use units such that a — 1 , fi = 1 , 
e = 1, and Boltzmann constant kp = 1. We will analyze 
bipolaron systems in thermodynamic limit (up to L = 32 
for rj = 0.1) measuring energies in units of t ~ 80meV. 
US] We fix U = 2.0t. HO] 

Method The temperature range relevant for transport 
properties is such that uiBuik < T, where the adiabatic 
approximation for bulk lattice degrees can be considered 
valid. [TT] Therefore, the dynamics of intermolecular bulk 
modes can be treated as classical. The electron dynam- 
ics is strongly influenced by the lattice oscillations which 
are described by the distribution function P ({x^}), which 
is self-consistcntly calculated through a Monte-Carlo ap- 
proach. [TT] At a fixed configuration of bulk displacements 
{xi}, Eq.@ shows that the electron hopping between 
nearest neighbor sites is not homogeneous and depend on 
the specific pair. 

At a fixed configuration {x;}, the resulting inhomogc- 
neous model can be accurately diagonalizcd using the fol- 
lowing basis states 

Ui({xi})\n qi ,.,n qL >\i,j> Sz =o, (8) 

where the unitary transformation U ({xi}) 



u ({ x i}) = ex P 
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(9) 

provides an el-ph treatment equivalent to a general- 
ized Lang-Firsov approach. [27] For long-range el-ph in- 
teractions, this approach is very accurate in the in- 
termediate to strong coupling regime. |28] The param- 
eters f q ({xi}) can be self-consistently calculated or, 
without losing accuracy, parametrized as f q ({x;}) = 
a({xi}) (M q /ui n ter)/(b({xi})cos(q) + 1), with a({xi}) 
and b ({xi}) configuration dependent variational quanti- 
ties. [33] We point out that, after the diagonalization at 
fixed displacements {xi}, the search for the best varia- 
tional parameters f q ({x;}) is made through the minimiza- 
tion of the free energy F. This procedure is based on inten- 
sive exact diagonalizations. Through the variational pa- 
rameters f q {{x{\), interface phonons are displaced from 
their equilibrium position to a distance proportional to 



the el-ph interaction. In Eq. 



> denotes the 



phononic basis in the occupation number representation, 
and \i,j >s 2= o= c lt c ]j.|0 >7 with |0 > vacuum state, is 
the basis state for the electron pair with total spin com- 
ponent S z = 0. We consider this Hilbcrt subspace since, 
for large el-ph coupling, the bipolaron is in a singlet state 
(S = S z = 0). 

In order to clarify the nature of the obtained solutions, 
some correlation functions such as the double occupation 
d =< J2i n i.t n i,i > an d the nearest neighbor quantity 
p =< J2i n i,t n i+i,i > are calculated. The main aim of 
this work is the mobility as a function of the temperature. 



The mobility \x ({x;}) is determined starting from the real 
part of the conductivity Re[a(oj)} ({x/}), taking the limit 
of small frequency and dividing for the particle density 
2/L of the system. In the linear regime, the real part of 
the conductivity is derived from the Kubo formula |30] 



Re[<r(w)] ({x,}) 
where 

3 = ' 



[1 - e-fr") 1 



2uj 



(10) 



(11) 



is the current operator, hj t s ({xi}) = t — aS(xj + s — Xj) 
denotes the generalized hopping and 5 = —1,1 indicates 
the nearest neighbors. 

If an observable O depends on the displacements {x^}, 
first, one makes the average (O) ({xj}) over the eigen- 
states and eigenvectors of the corresponding inhomogc- 
neous model, then over the distribution P({x,}) making 
the integral 
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(12) 



by means of the Monte-Carlo procedure. [UJ 

Results First, we analyze the stability of bipolarons com- 
paring the free energy of the two-particle problem i*2P 
with that of two single polarons Fp. The bipolaron is sta- 
ble if the difference A/ = F^p — 2Fp < 0. In the upper 
panel of Fig. 1, we report the free energy difference A/ for 
different values of interface and bulk el-ph couplings as a 
function of the temperature. By means of the correlation 
functions d and p, we can follow the crossover between 
large polaron (LP), large bipolaron (LB), and small bipo- 
laron (SB) states. The adjective large means extended 
range of lattice distortions, which implies that, in f q , a 
is smaller than unity and b is different from zero. The 
adjective small refers to local lattice distortions, and cor- 
responds to a = 1 and 6 = 0. In Fig. [TJ we show that 
the SB is characterized by the double occupation d = 1 
and the nearest neighbor correlation p = pointing out 
that the SB is made of two particles with spin up and spin 
down on the same site. On the other hand, the LB solu- 
tion exhibits a moderate value of d and p, both decreasing 
with temperature. As expected, the wave-function of LB 
spreads out over many sites, and the spreading gets larger 
with increasing temperature. [28] 

Surprisingly, bulk el-ph coupling tries to change the 
structure of the LB and SB solutions reducing their stabil- 
ity. As shown in the upper panel of Fig. [TJ the effect of the 
bulk coupling (Xp = 0.1) forbids the formation of bound 
states at Xj = 0.7. Apparently, the bulk coupling alone 
is not able to bind the two particles. The bulk SSH-like 
interaction reduces the ratio between the particle average 
kinetic energies with and without e-ph coupling, tending 
to enhance the localization of the particle and decreasing 
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Fig. 1: The free energy difference A/ (in units oii), the double 
occupation d and the nearest neighbor correlation function p 
as a function of the temperature (in units of t) for different 
values of interface (A/) and bulk (As) el-ph couplings. (LP) 
stands for large polaron, (LB) large bipolaron, and (SB) small 
bipolaron. 



the free energy Fp. With increasing the bulk el-ph inter- 
action, similarly, the electron pair reduces the free energy 
Fb- However, the bulk el-ph coupling tends to extend the 
spatial correlation of the local SB solution decreasing d 
and increasing p. On the other hand, in the LB regime, 
the bulk el-ph interaction tends to enhance both d and 
p. Therefore, the bulk interaction tends to decrease the 
pairing range for the LB going towards bond correlations. 

In Fig. [3 we show the phase diagram of the bipolaron 
for Xb = and Xb = 0.1 in the plane A/ vs. temper- 
ature. We notice that, in the limit of a single particle, 
the crossover (dotted line) from LP to small polaron (SP) 
occurs for high values of A/ in agreement with previous 
studies. [T()1[T7] We stress that the crossover line from LP 
to LB is shifted towards larger A/. The crossover line 
between LB and SP is less slightly shifted, but always to- 
wards larger interface couplings. 

The main result of this study is related to the analysis of 
the transport properties. As shown in Fig. [3[ the mobility 
has been calculated for different values of Xi at Xb =0.1 
focusing on the regimes of LP, LB, and SB. For Xi = 0.4 
(LP), the mobility fi scales as a power-law (/x ~ T~ 7 , with 
7 = 1.72). This behavior is due to the electron diffusion 
by bulk modes. QT] The main effect of the interface el-ph 
coupling is to reduce the exponent 7 starting from values 
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Fig. 2: Phase diagram of the bipolaron in the plane tempera- 
ture (in units of t) vs. interface el-ph coupling A/ for As = 
(dash line) and As =0.1 (solid line). LP stands for large po- 
laron, LB for large bipolaron, SB for small bipolaron, and SP 
for small polaron. The dotted line marks the crossover between 
large (LP) and small polaron (SP) for As = in the limit of a 
single particle. 



close to 7 ~ 2. For these regimes of parameters, the role 
of the interface coupling is similar to that of bulk intra- 
molecular local modes investigated in a recent paper. [32] 

For Xj = 0.9, we observe a crossover between two 
regimes. Below T = 0.25 (see Fig. [TJ, the LB is stable. 
We still find a power law (fj, ~ T^ 1 ) as in the case of LP, 
but with the reduced exponent 7 ~ 0.7. Above T = 0.25, 
the mobility still decreases with temperature, but the ex- 
ponent 7 is further reduced. Actually, the LP is on the 
verge of the crossover to SP which occurs at temperatures 
much larger than O.'St. Our findings for intermediate cou- 
pling regime are relevant for the experimental situation in 
which the gate is poorly polarizable. [4l-[6jl33j . 

The mobility in SB regime is shown in Fig. [3] for 
Xi = 1.2. It is strongly reduced and changes from a tiny 
metallic behavior at low temperatures to an insulating ac- 
tivated regime at high temperatures. Actually, at high T, 
p, goes as exp(— A/T), with A ~ 0.97t. Therefore, the 
bipolaron mobility naturally justifies activation gap of the 
order the hopping electron energy t for realistic values of 
the interface el-ph coupling. We stress that, in the limit of 
a single particle, [T7] for A/ = 1.3, A is of the order of 0.25i, 
and, for the unrealistic value A/ = 2.1, A is only about 
OAlt. Finally, we note that the order of magnitude of the 
mobility is comparable with experimental data in rubrene 
OFETs. [5] Actually, the mobility is expressed in natu- 
ral units, that is in terms of ^0 = ea 2 /h ~ 7cm 2 /(V ■ s), 
taking a = 7 A. [8]. Therefore, the mobility is about 0.1 
cm 2 /(V ■ s) at room temperature, that is two orders of 
magnitude less than that of OFETs exposed to vacuum. 

Conclusions and discussions In summary, OS gated 
with polarizable dielectrics provide conditions favorable 
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Fig. 3: The mobility fj, (in units of )j,o = ea 2 /K) as a function 
of the temperature (in units of t) for different strengths of the 
interface el-ph coupling: \i = 0.4 for the large polaron (LP), 
Aj = 0.9 for the large bipolaron (LB), \i — 1.2 for the small 
bipolaron (SP). The quantity 7 is the exponent of the power- 
law T -7 , while A is the activation energy for the hopping law 
exp(-A/T). 



for the el-ph induced bipolaron formation. We have exam- 
ined a model which includes both the bulk el-ph coupling 
and the long-range interaction of the carriers with the gate 
modes. The crossover to small bipolarons occurs for real- 
istic el-ph coupling strengths much smaller than those for 
the small polaron formation. However, within the small 
bipolaron regime, the mobility is characterized by an acti- 
vated behavior with gaps much larger than those obtained 
in the small polaron regime. The results have been dis- 
cussed in connection with rubrene field-effect transistors. 

In this work, we have assumed a simpler interface 
phonon spectrum and el-ph coupling, in comparison with 
models in high dimensions. |13j However, the interface el- 
ph coupling shows a long-range behavior similar to that in 
the more realistic two-dimensonal case. Actually, an im- 
portant role is played by interface el-ph coupling: in fact, 
Winter is of the order of or larger than 1, while Xsuik is 
of the order o/O.l. The features due to bulk coupling are 
limited to small changes in the bipolaron phase diagrams 
and the tiny coherent contribution to mobility at very low 
temperatures. 

This work is focused on the regime of low density since 
there are a few examples of highly doped OFET on solid 
substrates. [T1IT2"] With increasing density, in the regime of 
strong el-ph interaction, the SB phase is found to be stable 
upon clustering. [28j However, in the intermediate regime, 
the situation can change upon slightly increasing doping. 
Actually, a clustering of polarons or bipolarons could take 
place, and there should also be the possibility of a phase 
separation between very large undoped regions with large 
lattice deformations and tiny regions at substantial fill- 
ing with small deformations. Therefore, the increase of 



d isity involves peculiar many-body effects which require 
d ip investigations. 

11 addition to transport properties, the possibility 
tc distinguish between polaron and bipolaron formation 
c< ild come from optical spectroscopy and electron spin 
n onancc experiments. |28j However, there are not so 
n ,ny experiments with single crystal rubrene on highly 
p arizable dielectric gates. To our knowledge, the only 
n asurements of optical properties of rubrene have been 
p 'formed with gates made of the poorly polarizable 
p -ylene [3J] where polaron/bipolaron signatures do not 
ei erge. Similar situation occurs for experiments based on 
el ctron spin resonance. [2] Therefore, the bipolaron pro- 
p sal is expected to be more firmly confirmed by different 
p )bes not only in rubrene but also in other OS. 
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